Effective model of the electronic Griffiths phase 
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We present simple analytical arguments explaining the universal emergence of electronic Griffiths 
phases as precursors of disorder-driven metal-insulator transitions in correlated electronic systems. 
A simple effective model is constructed and solved within Dynamical Mean Field Theory. It is shown 
to capture all the qualitative and even quantitative aspects of such Griffiths phases. 
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I. INTRODUCTION 

The recent discovery of a number of heavy fermion 
materials with non-Fermi liquid (NFL) thermodynamic 
and transport properties has been followed by a signif- 
icant theoretical effort to understand the origin of the 
NFL behaviori In the cleaner systems, the proximity 
to a quantum critical point appears to dominate many 
of the observed propertiesiSiiMiSiSiL&iLiii In other heavy 
fermion systems disorder seems to play a more essen- 
tial role and appears to be crucial for understanding the 
NFL behavion 11 : 12 ! 13 : 14 ! 15 : 16 ! 17 Many experiments can 
be explained by the disordered Kondo model jii which 
has recently been put on a much stronger microscopical 
foundationi 1 ?i 1 ?i 2( i ) 

The emergence of electronic Griffiths phases in mod- 
els of correlated electrons has been establishedi&iiLSiiSi 
as a universal phenomenon, within a class of extended 
("statistical") Dynamical Mean Field Theory (DMFT) 
approaches. 21 This staiDMFT method provides an ex- 
act (numerical) treatment of localization in the absence 
of interactions, and reduces to the standard DMFT 
equations 2 ^ in the absence of disorder. When both in- 
teractions and localization are present, non- Fermi liq- 
uid behavior emerges universally^* as a precusor of a 
disorder-driven metal-insulator transition, due to a very 
broad distribution P(Tk) of local Kondo temperatures. 
This distribution has a Iow-Tr- tail of the form P(Tk) 
~ Ty* -1 , independent of the microscopic details or the 
specific form of disorder. The exponent a = a(W) is 
found to be a smooth, monotonically decreasing func- 
tion of the disorder strength W, and the NFL behavior 
emerges for W greater than the critical value W n fi cor- 
responding to a < I, when P(Tk) becomes singular at 
small Tk- As in other Griffiths phases, the thermody- 
namic and transport properties in this NFL region are 
dominated by rare events, which in this model correspond 
to sites with the lowest Kondo temperatures. 

In this paper, we show that the same behavior is found 
in a simpler, standard DMFT version of the model with 
a judicious choice of bare disorder. We should emphasize 
that localization is not present in this effective model, 
but the Griffiths phase emerges in qualitatively the same 



fashion as in the above more realistic calculations. We 
discuss how the specific disorder distribution which is 
hand-picked in the effective model is dynamically gener- 
ated by fluctuation effects within the statDMFT formu- 
lation, elucidating the origin of the universality of the 
Griffiths phase behavior. In addition, the simplicity of 
this DMFT effective model makes it possible to describe 
all the qualitative features of the solution using simple 
analytical arguments, thus eliminating the need for large 
scale numerical computations in the description of the 
electronic Griffiths phase. This may be crucial in order 
to address more complicated issues, such as the role of 
additional Ruderman-Kittel-Kasuya-Yosida (RKKY) in- 
teractions in disordered Kondo alloys^ 

This paper is organized as follows. Sec. ITT1 introduces 
the effective model for the electronic Griffiths phase as 
a DMFT model with a Gaussian distribution of random 
site energies. This model is solved analytically in the 
Kondo limit in Sec. IIIII and numerically in Sec. IIVI The 
arguments explaining the universal aspects of the form 
of the renormalized disorder are presented in Sec. Ivl 
Sec. IVII establishes a connection with the Griffiths phase 
in a single band Hubbard model, and Sec. IVIII contains 
our conclusions. 



II. MODEL 

We consider the Anderson lattice model where the dis- 
order is introduced by random site energies e, in the con- 
duction band, as given by the Hamiltonian 
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where fj a and Cj a are annihilation operators for /- and 
conduction electrons, respectively. V is the hybridization 
parameter, and Ef is the /-electron energy. We assume 
U — > oo, and choose a Gaussian distribution of random 
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site energies for the conduction band 

P( Si ) = (27TW 2 )- 1 / 2 exp{-~4/W 2 }. (2.2) 

In Sec. E] we will explain how this particular disorder 
distribution comes out naturally from the more generic 
statDMFT approach. 

To solve these equations, we use the DMFT 
approach^ which is formally exact in the limit of large 
coordination. We concentrate on a generic unit cell j, 
containing an /-site and its adjoining conduction elec- 
tron Wannier state. After integrating out the conduc- 
tion electron degrees of freedom, we obtain the effective 
action for the /-electron on site j 

S mp (j) = £ f dr f dT'fl(T)[5{T-T'){d T +E f ) 
a Jo Jo 

+ A fj (T-r')}f j(7 (T'). (2.3) 

Here, the restriction of no double /-site occupancy is 
implied. The hybridization function Afj between the /- 
electron and the conduction bath A c is given by 



A fj (iu n ) 



V 2 



Ej + fj,- A c (iu) n ) 



(2.4) 



The self-consistency condition for the conduction bath 
(cavity field) assumes a simpler form for the semi-circular 
model density of states^, which we use for simplicity. All 
the qualitative features of our solution are independent 
of the the form of the lattice, and the quantitative results 
depend only weakly on the details of the electronic band 
structure. For this model A c (iu> n ) = t 2 G c {iuj n ), where 
G c {iuj n ) is the disorder-averaged Green's function of the 
conduction electrons, and the self-consistency is enforced 
by 

G c (iUJ n ) = ([iaJn-Ej+li-^GciiLUr^-^j^aJn)]" 1 ), (2.5) 



where 



V 2 



iui n - E f - T,fj(iu! n ) ' 



(2.6) 



and Tifj is the /-electron self-energy derived from the 
impurity action of Eq. I|2.3[) . From a technical point of 
view, within DMFT the solution of the disordered An- 
derson lattice problem reduces to solving an ensemble 
of a single impurity problems supplemented by a self- 
consistency condition. 

We will solve the system of Eqs. (|2.3[) - (I2.6|) at zero 
temperature using the slave boson mean field theory 
approachi22i2£ This approximation is knowni£ii2*2£ to re- 
produce all the qualitative and even most of the accu- 
rately quantitative features of the exact DMFT solution 
at T = 0. It introduces renormalization factors (quasi- 
particle weights) Zj and renormalized /-energy levels e fj , 
which are site-dependent quantities in the case of a dis- 
ordered lattice. These parameters are determined by the 



saddle-point slave boson equations (see Ref. [2(| for more 
details) which, on the real frequency axis, assume the 
form 
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(2.7) 



Eq. H2.6JI in this case becomes 



(2.8) 
(2.9) 



III. ANALYTICAL SOLUTION IN THE KONDO 
LIMIT 



Before presenting a numerical solution of the slave bo- 
son Eqs. (|2.7() - (|2.8|l supplemented by the self-consistency 
condition of Eq. I|2.5|) , we will solve these equations ana- 
lytically in the Kondo limit for a given conduction bath. 
A comparison with the numerical solution will show that 
the self-consistency does not qualitatively change the an- 
alytical results. 

The slave boson equations simplify in the Kondo limit 
Zj — > 0. The integral in Eq. (|2.7I) is dominated by the 
low-frequency region, and the frequency dependence in 
A c and Afj can be neglected. Therefore, after integra- 
tion 



-ZjReiAfM] 



(3.1) 



where, for simplicity, we took a semi-circle conduction 
bath with /i = 0. In the integral of Eq. 12.8fl . the fre- 
quency dependence of Afj can also be neglected. Intro- 
ducing the energy cutoff D and using Eq. (|3.1|l we obtain 



e 2 + (7Tt 2 Po ) 2 
D- 



Zo- 



nt 2 Po V 2 
) + {*t 2 Po f 

(7Tt 2 Po ) 2 



-7T 2 i 2 p l .J -e) ft 2 p J Ej /2t 2 Po 



■ exp 



t 2 Po J 



2e, 



(3.2) 



Here, p is the density of states (DOS) of the conduction 
electrons at the Fermi level, J = 2V 2 /\Ef\, and Z = 
Z (ej = 0) . The Kondo temperature is proportional to 
the quasi-particle weight, Txj — irV 2 p Zj. In the limit 
Ej ^ J/2 and neglecting a weak site-energy dependence 
in the pre-factor, we obtain 



Tj 



T -1/Xj 



where the site dependent coupling constant is 

t 2 p J 
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(3.3) 



(3.4) 
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and TS- is the Kondo temperature in the clean limit (for 
£j = 0). From these equations, we can immediately 
find the desired distribution of local Kondo temperatures 
P(T K ) = P(e(T K ))\de/dT K \, which (up to a negligible 
logarithmic correction) is given asymptotically by 



P{T K ) cx (T K /T° K Y 



with 



a(W) 



t 2 p J 
2W 2 ' 



(3.5) 



(3.6) 



This expression is one of the central results of this pa- 
per. It has exactly the form expected for a Griffiths 
phase, where the exponent characterizing the local en- 
ergy scale distribution assumes a parameter-dependent 
(tunable) form. 

To show how the NFL behavior appears due to the 
singularity in P(Tk), we use the standard expression due 
to Wilson for the magnetic susceptibility 27 



X(T,T K ) = 



C 



T + oTk ' 



(3.7) 



which is an excellent approximation for a single Kondo 
impurity. Here C and a are constants. In the disor- 
dered case, we can split the average susceptibility x(T) = 
P{T K )x{T,T K )dT K in a regular "bulk" part 



rOO 

Xr (T)= / P(T K ) X (T,T K )dT K , 

J A 

and a potentially singular part 
Xs(T) = C x 



A rpa—1 

T + aT K 



dT K , 



(3.8) 



(3.9) 



coming from the tail with low Kondo temperatures (A 
is a crossover scale). At weak disorder, the exponent a 
is large and the distribution P(Tr-) is regular, x(0) = 
Xo + C2/{a — 1), but NFL behavior emerges once a ^ 1, 
which corresponds to 



W > W nf l = y/t 2 Po J/2. 



(3.10) 



For a = 1 the magnetic susceptibility has a logarithmic 
divergence, x(T) oc ln(l/T), characteristic of marginal 
Fermi liquid behaviorj2& while for a < 1 a power law 
divergence is obtained, x(T) cx T Q_1 as T — > 0. The 
same singularity also leads to an anomalous behavior in 
the transport properties, as shown in detail in Refs. 1261 
and El 



IV. NUMERICAL RESULTS 

In the above derivation we ignored the fact that the 
conduction bath A c has to be self-consistently deter- 
mined. This will also produce particle-hole asymmetry 




0.0001 



FIG. 1: Distribution of the local Kondo temperatures. The 
disorder ranges from moderate W = 0.1 to strong W = 0.4. 
We used V = 0.5, E f = -1, and /i = -0.1. 



and an asymmetric distribution of Kondo temperatures 
Tjcj. A nonzero chemical potential will further increase 
this asymmetry. However, the numerical solution we ob- 
tained using the slave boson approximation at zero tem- 
perature shows that the essential physics described by 
Eqs. I|3.3() - (|3.6I) remains qualitatively correct. The dis- 
tribution of local Kondo temperatures in the asymptotic 
limit is indeed a power law, P(Tr-) ~ Tx Q_1 , where the 
exponent a is a decreasing function of disorder. 

Fig. shows the distribution P(Tr-) for several values 
of the disorder distribution strength W. For the param- 
eters that we here use, the system is close to the Kondo 
limit, and the Kondo gap of the clean system is approx- 
imately 0.04 (in energy units of the half bandwidth of 
bare DOS). The NFL behavior appears for W > 0.14. 
We note that in the NFL regime the power law behavior 
appears already for the site energies Ej which deviate only 
moderately from the mean (zero) value. In other words, 
the asymptotic behavior is established well before we at- 
tain very rare realizations of Sj which belong to the tail 
of the Gaussian distribution. For example, for W = 0.3, 
sites with Ej > 0.4 (which correspond to Z < 0.01) are 
already in the power-law regime. 

According to the simplified derivation from Scc lIIII the 
exponent a is inversely proportional to W 2 . The numeri- 
cal results shown in Fig.[21confirm such behavior for weak 
and moderate disorder. For strong disorder there appear 
some deviations from this formula, which can be ascribed 
to the dependence of the DOS at the Fermi level on the 
disorder strength. 

Before we present arguments which justify our effec- 
tive DMFT model approach, let us make a direct com- 
parison with the statDMFT results from Ref. |2(]. In 
this approach, very broad distributions of local Kondo 
temperatures are generated for arbitrary distributions of 
bare disorder. In particular, even if the bare distribution 
is bounded, sites with arbitrarily small Kondo tempera- 
tures will exist, and their distribution will have a power 
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FIG. 2: Inverse power law parameter a 1 as a function of W 2 . 
For weak and moderate disorder this dependence is linear. 
Here V = 0.5, E f = -1, and fi = -0.2. 
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FIG. 3: Power law exponent a as a function of disorder 
strength measured by the standard deviation W . Full lines 
are the effective model results, and dashed lines are the 
stafDMFT results. The hybridization V is taken to be 0.5 
and Ef = — 1. The chemical potential is fi = —0.1 (upper 
lines), and —0.5 (lower lines). The onset of NFL behavior 
occurs at a = 1. 



law tail. This is a consequence of the spatial fluctuations 
of the conduction electron cavity field, as we discuss in 
detail in the next Section. In Fig. [3] we compare the val- 
ues of the exponent a for the effective model with Gaus- 
sian disorder of variance W 2 , and the statDMFT results 
obtained for a bounded uniform distribution of bare dis- 
order with the same variance. Remarkably, not only does 
the electronic Griffiths phase emerge in the same fashion, 
but the numerical values of disorder strength determin- 
ing the onset of NFL behavior are also almost the same. 
The comparison is made for two different values of the 
chemical potential. As we move further away from half- 
filling by changing the chemical potential, the critical 
value W n fi decreases. That is expected since a should 
be proportional to the bare (noninteracting) DOS at the 
Fermi level. 



V. THE ROLE OF SPATIAL FLUCTUATIONS 
AND THE FORM OF RENORMALIZED 
DISORDER 

In this Section we explain the universal aspects of the 
emergence of the electronic Griffiths phase within the 
more generic statistical DMFT. In particular, we show 
how the Gaussian tails in the distribution of renormal- 
ized disorder appear for an arbitrary form of the bare 
disorder. Moreover, we present arguments showing that 
the Griffiths phase appears generically as a precursor of 
the Mott- Anderson metal-insulator transition. 



A. Universality of the renormalized disorder 
distribution 

In the above DMFT formulation, we had to choose 
a special form of disorder distribution in order to ob- 
tain the desired power-law distribution of Kondo tem- 
peratures. Had we chosen a different distribution, the 
results would not have held. For example, for a bounded 
distribution of site energies, there would always be a 
minimum value of the Kondo temperature, and thus no 
power-law tail. On the other hand, from numerical sim- 
ulations of lattices with finite coordination, it has been 
established that the emergence of the Griffiths phase is a 
universal phenomenon^ Why? To understand the rea- 
son for this, we note that for finite coordination (as op- 
posed to the DMFT limit) the cavity bath A c is not 
self-averaging, but is a site-dependent, random quantity 
A C j . In this statDMFT formulation, the local conduction 
electron Green's function is given by 

G c j(iu) n ) = — r ——, r, (5.1) 

lLU n - Ej + fj, - A cj {lLU n ) - <Bj[lW n ) 

where $j describes the local scattering of the conduction 
electrons off the /-shell at site j, and is given as before 
by Eq. 

For weak disorder, the corresponding fluctuations are 
small, and we can separate 

A cj = A a c v + SA cj . (5.2) 

In the following, we compute the distribution function for 
the fluctuations of the cavity field, which will lead to the 
renormalized form of the disorder distribution function. 
The renormalized site energies can be defined as 

e j =e j +5A^ j , (5.3) 

where 8A^ — Re [SA c j(w — 0)]. We stress that the cav- 
ity fluctuations are present for a general finite coordina- 
tion electronic system in the presence of disorder of any 
kind. In particular, the disorder in hybridization param- 
eters Vj, or local /-energy levels Efj, will induce fluctua- 
tions in the local DOS even if random site energies Ej in 
the conduction band are absent. Furthermore, as we ar- 
gue in the next subsection, the renormalized distribution 
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P(£j) will have universal Gaussian tails even if the bare 
distribution P(sj) is bounded. Note that 6A C j has a real 
as well as an imaginary part SA 1 ^, due to the fact that 
fluctuations locally violate particle-hole symmetry. How- 
ever, we show in the Appendix that S A^- fluctuations, at 
least when treated to leading order, do not produce sin- 
gular behavior in P{Tk) and therefore can be neglected 
when examining the emergence of the electronic Griffiths 
phase. 



B. The Gaussian nature of the renormalized 
distribution 

From detailed numerical studies it has become clear 
that the onset of the Griffiths phase in disordered An- 
derson lattices generally occurs already for a relatively 
moderate amount of disorder. In this limit, the relevant 
distributions are determined essentially by the central 
limit theorem, therefore resulting in a Gaussian form of 
the tails for P(sj)- This is precisely what is needed to 
justify the DMFT effective model, where such Gaussian 
tails are assumed from the outset. 

Before engaging in more precise computations of these 
distributions, it is worth pausing to comment on the 
physical validity of the assumed Gaussian statistics, i.e. 
the relevance of the central limit theorem in the cases 
of interest. Quite generally, if a certain quantity can be 
represented as a sum of a large number of independent 
random variables, then the central limit theorem tells 
us that the resulting distribution will be Gaussian, ir- 
respective of the specific form of the distributions of the 
individual terms in the sum. In our case, the fluctuations 
of the local cavity field result from Friedel oscillations of 
the electronic wave functions, induced by other impuri- 
ties which may lie at a relatively long distance from the 
given site. This is a result of the slow (~ R~ d ) decay of 
the amplitude of the Friedel oscillations in d dimensions, 
where R is the distance from the impurity site. The sit- 
uation is very reminiscent of the Weiss molecular field 
of an itinerant magnet, where the RKKY spin-spin in- 
teractions have a long range character for the very same 
reason, being as they are a reflection of similar Friedel 
oscillations. Furthermore, as we will explicitly show, the 
leading corrections (to order O (W 2 )) at weak disorder 
take the form of a linear superposition of contributions 
from single impurity scatterers, and thus of a sum of 
independent random numbers, for which we expect the 
central limit theorem to hold. 

To obtain the precise form of this distribution, it there- 
fore suffices to compute the variance 



function in a random medium. At zero frequency for this 
site 



al = ((5A*) 2 > , 



(5.4) 



to leading order in disorder strength. To compute the 
fluctuations 6A C j at weak disorder, we note that the cav- 
ity field A c j can be computed if we consider a particular 
site with Ej = (call it site 0), and compute its Green's 



A co = fx - 1/G 00 . 
The corresponding variation is 

SA co = SG 00 /(G 00 ) 2 . 



(5.5) 



(5.6) 



We still need to compute the fluctuation SG 00 , which can 
be expanded in powers of the random potential Sj . In do- 
ing this, we have ignored the interaction renormalizations 
of the random potential for conduction electrons. We will 
return to re-examine this effect in the Appendix. Note, 
however, that in the absence of interactions in the envi- 
ronment of a given site, the following expressions provide 
the exact leading contributions at weak disorder. 
To leading order, we can write 



This gives 



where 



6G 00 = J2 s j(Goj) 2 + 0(e 2 ). 



o% = CW 2 + 0{W 4 ), 



Re( G -) 2 



(Goo)' 



(5.7) 



(5.8) 



(5.9) 



The Green's function sum C will numerically depend on 
the lattice geometry, but will generally be a dimensionless 
number of order one. 

The distribution of renormalized disorder e = e + SA^ 
is given by a convolution of the distributions Pi(s) and 



P(e) 



du P x (e - u)P 2 {w), 



(5.10) 



If the bare distribution is bounded, (e.g. a standard 
"box" distribution), then Gaussian tails will emerge due 
to the fluctuations in SA^, and the "size" of the tails will 
be determined by an effective disorder corresponding to 
Wj// = WVC. Here, the superscript indicates that 
we ignored the interaction renormalizations. In the Ap- 
pendix, we argue that the effective scattering potentials 
$j will further renormalize the disorder distribution, but 
the Gaussian tails will remain as its generic feature. 

Now we present numerical results which provide strong 
evidence for the universality of the renormalized disorder. 
Fig.EJshows the results obtained within the statDMFT 18 
for uniform and binary disorder distributions with the 
same standard deviation W = 0.1. As anticipated by 
Eq. I|5.8I) , the fluctuations of the cavity bath acquire an 
approximately Gaussian form with the same variance for 
both bare disorder distributions, panel (b). The renor- 
malized disorder distribution P(e) exhibits long tails, 



6 



1(a)' A ' A 


■ i i 

uniform _ 

binary . 


i 1 i 


1 1 
1 1 
1 1 


1 i 


1 1 
I I I 


- 

i 1 



-0.5 



0.5 



1 1 




1 1 1 


1(b) 


uniform - 






i binary " 




I 





-0.5 



0.5 



SA' 



. (c) 


1 1 


1 1 1 1 






^"Sv. uniform _ 






binary 


i r"^"^ i 


• 





FIG. 4: StatDMFT results for the disorder distributions. Uni- 
versal Gaussian-like tails appear in the cavity field fluctua- 
tions, panel (b), and the renormalized disorder distributions, 
panel (c), even though the bare disorder distributions are 
bounded (uniform and binary), panel (a). We used V = 0.5, 
Ef = —1, and p — —0.5. 



panel (c), although the bare distributions are bounded, 
panel (a). These Gaussian-like tails are the main univer- 
sal feature of the renormalized disorder, and they are cru- 
cial for the appearance of the singular behavior in P(Tk) 
which leads to the formation of the Griffiths phase. 

StatDMFT results in Fig. [S] provide further illustra- 
tion of the universality. The upper panel shows that the 
distributions of Kondo temperatures for uniform and bi- 
nary bare disorder distributions with the same standard 
deviation W are qualitatively the same. The exponent 
a, which determines the slope of the distribution tails, is 
shown at the lower panel as a function of W . It depends 
very weakly on the particular form of disorder distribu- 
tion. 



C. Localization effects 

In the strict DMFT formulation a w t 2 p Q J/2W 2 , 
where p = p av is simply the (algebraic) average DOS of 
the conduction electrons, which therefore remains finite 
even at the localization transition^ If J is chosen to be 
large enough, the above seems to suggest that the Grif- 
fiths phase may not emerge before the electrons localize 
at W = W c ~ 1/ pav However, in a theory that includes 
localization, the Kondo spins do not see the average, but 
rather the typical DOS of the conduction electrons 
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FIG. 5: StatDMFT results for the distribution of Kondo tem- 
peratures for uniform and binary bare disorder distributions 
(upper panel). The lower panel shows the exponent a as a 
function of disorder strength. The results are qualitatively 
the same regardless of the specific form of the disorder distri- 
bution. We used V = 0.5, Ef — — 1, and p — —0.5. 



Thus, in the NFL criterion, Eq. (|3.10() . one should ac- 
tually replace p — > pt yp , a quantity that becomes very 
small (and eventually vanishes) as the Anderson transi- 
tion is approached, viz. 



p typ = A{W C - Wf, 
where A and (3 are constants. We thus get 

Wl fl = \At 2 J{W c -W nfi y. 



(5.11) 



(5.12) 



This transcendental equation cannot be solved in closed 
form in general, but an approximate solution can be 
found for W*-P/At 2 J < 1. In this case, the quantity 
8W = 1 — W n fi/W c will be small, and to leading order 
in SW 

W nfl =W C - (At 2 Jffl-WwW < W c . (5.13) 

Therefore, the Griffiths phase emerges strictly before the 
transition is reached. 
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VI. ELECTRONIC GRIFFITHS PHASE IN THE 
VICINITY OF THE METAL-INSULATOR 
TRANSITION 

Previous workSi has shown that the electronic Griffiths 
phase appears also in a single band Hubbard model, as 
a precursor to the Mott- Anderson metal-insulator tran- 
sition (MIT). Since the Hubbard model at half-filling is 
equivalent to the charge-transfer model^, of the MIT, we 
examine in this Section the appearance of the Griffiths 
phase within this model, which can be considered a ver- 
sion of the Anderson lattice model that we examined in 
our approach. 

The charge transfer model has been used to describe 
the Mott metal-insulator transition for various systems, 
including copper oxides »2£ It consists of a two band 
model, where one of the bands is narrow, and has large 
on-site interaction U (Copper d-band), while the other 
band is broad enough that electron-electron interactions 
can be neglected (Oxygen p-band). A disordered version 
of this model is also appropriate to describe the Mott- 
Anderson transitional in doped semiconductors such as 
Si:P. Here, the narrow band corresponds— to the impu- 
rity band of the Phosphorus donors, while the broad one 
is the conduction band of Silicon. This model is given ex- 
actly by the Hamiltonian of Eq. (|2.1fl , but supplemented 
by the constraint 



(6.1) 



which can be enforced by adjusting the value of the chem- 
ical potential. Here n fj and n c j are the average number 
of conduction and /-electrons on site j, and the overbar 
denotes the average over disorder. In the mean field slave 
boson approach, the average occupancy of the /-site is 
equal to 



1 



(6.2) 



As the /-electron energy level is decreased (\Ef\ in- 
creased), the occupancy of the /-sites becomes larger: the 
charge is "transferred" from the conduction band. The 
transition to the Mott insulator is found for sufficiently 
large \Ef\. At least within DMFT, this metal-insulator 
transition has the same character as the more familiar 
Mott transition in a single band Hubbard model. As an 
illustration, we show on Fig.[S]the number of conduction 
electrons per site, n c = Z, as a function of Ef, in the 
clean limit. 

We have solved our effective model in the parameter 
regime relevant to the approach to the Mott-Anderson 
transition, and the results demonstrate the emergence of 
an electronic Griffiths phase in the same fashion as for 
the Anderson lattice model, consistent with statDMFT 
results Here, the /-level energy in Fig. [7J is mea- 
sured with respect to the middle of the conduction band, 
and not with respect to the chemical potential as in 
Sec. IIVI For the parameters used in Fig. [7J the system 
is in the mixed valence regime, not in the Kondo limit, 
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FIG. 6: Number of conduction electrons per site as a function 
of the /-level energy Ef. The approach to the insulating 
phase is linear, n c ~ Ef — Ef, where Ef is the critical value 
for the /-level energy. The hybridization V = 0.5, and Ef is 
measured with respect to the middle of the conduction band. 




FIG. 7: Distribution of local Kondo temperatures for several 
levels of disorder. The total number of electrons per unit cell 
is fixed to 1. The disorder ranges from W = 0.2 to 0.38. 
We used V = 0.5 and Ef — —1.3, where Ef is measured 
with respect to the middle of the conduction band. The inset 
shows the exponent a as a function of the disorder strength. 
The NFL phase occurs for a < 1. 



and stronger disorder is needed for the appearance of 
the NFL electronic Griffiths phase, again in close agree- 
ment with statDMFT resultsiSi These results demon- 
strate that our effective model proves capable to describe 
the emergence of the electronic Griffiths phase as a uni- 
versal phenomenon in correlated electronic systems with 
disorder. 



VII. SUMMARY AND OUTLOOK 

In this paper we have identified an analytically solvable 
infinite range model, which captures the emergence of 
electronic Griffiths phases found within the more generic 
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statDMFT approachesi 18 i 20 i 21 In this effective model, a 
specific distribution of disorder is postulated, leading to a 
power-law distribution of local Kondo temperatures and 
NFL behavior for sufficiently strong disorder. We have 
also presented arguments explaining how this specific 
form of randomness is universally generated by renor- 
malizations due to disorder-induced fluctuations of the 
conduction bath. In this sense our effective model can 
be regarded as a (stable) fixed point theory of electronic 
Griffiths phases. 

The main motivation for introducing this effective 
model lies in its simplicity, allowing an analytical solu- 
tion, and thus providing further insight into the mecha- 
nism for the emergence of the electronic Griffiths phase. 
Nevertheless, an essential ingredient is still missing from 
our Griffiths phase theory, namely the inter-site RKKY 
interactions between Kondo spins. According to the 
existing picture, all the spins with Tk < T will not 
be Kondo screened, thus providing a large contribu- 
tion to thermodynamic response. These decoupled spins 
will, however, not act as free local moments, but will 
feature dynamics dominated by frustrating inter-site 
RKKY interactions. Recent experiments on disordered 
Kondo alloys indeed seem to suggest the presence of low 
temperature glassy dynamics with a negligible freezing 
temperature^*^ The simplifications introduced by our 
effective model open an attractive avenue to incorporate 
both the Kondo effect and the RKKY interaction in a 
single theory. This fascinating direction remains a chal- 
lenge for future work. 
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and we get 



where 



a] = CiW 2 + 0(W 4 ), 



Im 



-i 2 



{GqoY 



(A2) 



(A3) 



0*0) 



In this approximation, the quantity ^A^- has a Gaussian 
distribution, and we find 



P[T K ) ~ T^ 1 exp 



J 2 
'27rV 2 



\n 2 (D/T K ) 



(A4) 



As we can see, because the "log" in the exponent has an 
extra power of two, this distribution is log-normal and 
not power-law. Therefore the ^A^- fluctuations, at least 
when treated on the Gaussian level as we have done, do 
not lead to a singular P{Tk) distribution. Thus, to lead- 
ing order we can ignore these fluctuations when examin- 
ing the emergence of the electronic Griffiths phase. 



APPENDIX B: INTERACTION 
RENORMALIZATIONS 

In these estimates, we have omitted an important in- 
gredient, the fact that Kondo disorder itself will pro- 
duce additional scattering, i.e. disorder in the conduc- 
tion channel, which needs to be self-consistently deter- 
mined. As we have shown in previous work^S this re- 
sults in a distribution of effective scattering potentials 
$j, corresponding to the Kondo spins (note that in the 
uniform case, the $j-s are the same on all sites, result- 
ing in no scattering, but contributing to the formation of 
the Kondo gap). The resulting scattering, in the weak 
disorder limit again can be considered as a Gaussian dis- 
tributed potential of width 



APPENDIX A: FLUCTUATIONS IN SA' C] 

In Sec. V we have ignored the fluctuations in the imagi- 
nary part of the cavity function SA 1 ^. The corresponding 
contribution to the Iow-Tr- tail is sub-leading, as we now 
show. We need to focus on rare events that produce ex- 
ceptionally small values of the local conduction electron 
DOS p-j = -iG ci (0). Using Eq. ifOJl . and ignoring the 
fluctuations in <5A^, we see that low values for pj cor- 
respond to exceptionally high values for SA 1 ^. We there- 
fore need to compute the form of the high-JA^- tail of 
P(SA I c j). Just as for the real part, we can estimate the 
fluctuations of SA 1 ^ by calculating the second moment, 



<(^) 2 ) 



(Al) 



(Bl) 



Note however that this additional "Kondo" scattering 
does not enter directly (at site 0) in the solution of the 
local Kondo problem, since the local f-site "sees" the cor- 
responding c-site with the f-site removed. However, the 
presence of $j-s on all other sites (j ^ 0) does modify 
the form of SA^ which, therefore, has to be computed 
by including this additional scattering. At weak disorder, 
we expect 



<*?> 



(B2) 



where the constant C\ measures the response of the 
Kondo spins to the hybridization disorder. Note that 
W e ff enters here, since the <F,-s are obtained from the 
solution of local Kondo problems, which are determined 
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by the strength of the renormalized site disorder, as mod- 
ified by hybridization fluctuations. We therefore need to 
compute W e f / self-consistently, and we get 

W? ff = W 2 + C(W 2 + C 2 W 2 ff ), (B3) 

or 

Wef, = f^rW (B4) 

This reasoning, valid for weak bare disorder illustrates 
how the effective disorder is generated in the conduction 
band even if it originally was not there, or is enhanced 



due to additional Kondo scattering, if already present. In 
addition, these arguments illustrate how Gaussian tails 
are generated for the renormalized disorder, even if they 
are not introduced in the bare model. Of course, non- 
linear effects at stronger disorder cannot be accounted 
for in this simple fashion, which is especially true for 
the consideration of the additional scattering introduced 
by disordered Kondo spins. Nevertheless, the simple ar- 
guments that we presented illustrate how universality is 
produced by renormalizations due to cavity field fluctu- 
ations, and seem to capture the essential features of the 
emergence of the electronic Griffiths phase. 
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